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ABSTRACT 

The research program involved the modification and use of an existing pulse 
detonation engine (PDE) to investigate the detonability of a JP-10/air aerosol. The 
detonation of a JP-10 aerosol in air proved more difficult than was originally anticipated. 

The use of a small JP-10/oxygen pre-detonator to provide direct initiation results 
in a transition region with a geometric diffraction and compositional discontinuity. 
Propagation of a detonation into such a region is very complex but critical to the re- 
establishment of the detonation wave in the JP-10/air mixture. A high-speed camera was 
used to image the wave in the transition region and provide spatial information. High 
frequency pressure transducers were used along the combustor axis to determine wave 
speed. The ultimate goal was to determine the conditions required to ensure reliable re- 
establishment of a detonation wave in the JP-10/air aerosol mixture. 

Unfortunately, the confined planar JP-10/oxygen detonations in the pre-detonator 
were unable to transition into unconfined spherical detonation fronts in the JP-10/air 
aerosol. Furthermore, the ratio of main combustor diameter to pre-detonator diameter 


was too large to allow re-initiation of detonation at the main combustor wall. 
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I. INTRODUCTION 

A. PURPOSE 

This investigation was conducted as part of the Office of Naval Research program 
for fundamental research into pulse detonation engines (PDEs). During the past year, our 
goal was to determine the effects of fuel distribution on detonation wave stability and 
velocity in a JP-10/air aerosol. Testing apparatus consisted of a 1.5625-inch inner- 
diameter tube, or pre-detonator, surrounded by four similar fuel/air tubes, all of which 
opened into a 5.25-inch inner-diameter tube, or main combustor. (See Chapter IV and 
the Appendix.) A JP-10/oxygen aerosol was detonated in the pre-detonator for the 
purpose of initiating detonation in the JP-10/air aerosol in the main combustor. A 
parallel goal was to minimize the oxygen requirements of the pre-detonator. Previous 
work at the Naval Postgraduate School included detonation studies of ethylene/air, but 
for volume-limited systems (e.g., tactical missiles) liquid fuels are required for a good 
density impulse (pI,p) as well as avoiding the large tank mass required to contain 
pressurized gases. The ultimate objective of the program is to determine the viability and 
design requirements of a liquid-fueled PDE for tactical missile propulsion. 
B. BACKGROUND 

1. PDE Operation 

Detonations have been studied for a variety of purposes since the 19" century. 
The history will not be presented here but a very good summary is given by Bauer, 
Dabora, and Manson [Ref. 1, 2]. In a PDE, the detonation would be used to generate 


high-pressure combustion products. Figure 1 illustrates the pressure profile. 
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Figure 1. Pressure Profile within a PDE from Ref. [3] 


to determine the thrust. 


The pressure against the head wall (P3) of the PDE combustor can be integrated over time 
The operating cycle of a single-tube PDE is presented in Figure 2. 
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Figure 2. PDE Operating Cycle from Ref. [3] 


Initially, the fuel/oxidizer mixture is injected into the combustion chamber. As 
soon as the valves are shut, detonation is initiated at the head end. Since the detonation 
products flow in the direction of the wave, and the head end boundary is closed, 
rarefaction waves are generated. These reduce the pressure at the head wall (P3) below 


that immediately behind the detonation (P2). Head wall pressure can be calculated using: 
ay 
P= zat) P, (1) 
2y 

where y is the ratio of specific heats of the burned gases and P? is the pressure for the 
Chapman-Jouguet detonation condition (described in Chapter IT) [Ref. 4]. Once the 
detonation wave leaves the end of the tube, rarefaction waves are generated which travel 
towards the head end, reducing combustor pressure to ambient pressure. This serves to 
exhaust most of the combustion products. The cycle is then repeated. [Ref. 3] 

As the operating frequency increases, the pressure on the head wall will reach a 
quasi-steady value. Several factors govern the maximum cycle frequency including, 
physical tube size, fuel and oxidizer flow rates, valve response times, control program 
response times, and the time required for the combustion wave to transition to a 
detonation. Very high detonation repeatability is also required. Since deflagration is a 
much slower process than detonation, a failed detonation attempt could leave a slow 
combustion event still occurring in the tube when the fuel/oxidizer mixture 1s added to 
begin the next cycle. This new mixture could then begin burning before the tube is 


ignited. Thus a single detonation failure at a high operating frequency could prevent 


successful operation for several subsequent cycles. Multiple combustion tubes can also 
be used to increase the effective cycle frequency. 

As implied above, PDEs themselves are very simple. Aside from the valving, 
there are no moving parts, which should make PDEs inexpensive and reliable. Although 
not as simple as a ramet, the PDE has the advantage of being able to generate thrust from 
low subsonic through supersonic flight conditions. In addition, the cyclic nature of PDE 
operation makes it capable of near-instantaneous throttling over a wide thrust range. 

2. Initiation of Detonation 

There are two modes of detonation initiation: thermal initiation (or self-ignition) 
and direct initiation. In thermal initiation, a tube containing an explosive mixture is 
ignited at the closed end. The resulting combustion products have a specific volume 5-15 
times that of the reactants [Ref. 5]. This generates a series of compression waves, each 
one of which heats the unburned mixture to a slightly higher temperature. These higher 
temperatures lead to progressively higher flame speeds, which create stronger 
compression waves. In addition, sonic velocity in the unburned mixture is slightly higher 
for each successive compression wave, thereby allowing each to catch up with the initial 
one [Ref. 5]. Eventually the unburned gases become turbulent, accelerating the flame 
front even more until the compression waves coalesce into a shock of sufficient strength 
to ignite the unburned mixture [Ref. 5]. This reaction then supplies the energy required 
to sustain the shock, thus forming a detonation. This is the mode by which detonation is 


initiated in the pre-detonator. 


The other method of initiating detonation is direct initiation, in which a strong 
shock is directly applied to the explosive mixture. Due to the coupling between chemical 
and gas dynamic processes, the sensitivity of the reactants will dictate the required 
strength of the shock [Ref. 5]. The pre-detonator is used to provide the strong shock for 
direct initiation of detonation in the main combustor. This arrangement was chosen since 
the deflagration-to-detonation transition distance for a JP-10/oxygen aerosol is 
approximately 6-8 inches under optimum conditions [Ref. 6], while for a JP-10/air 


aerosol it would be orders of magnitude longer [Ref. 5]. 
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Il. DETONATION THEORY 
A. DETONATIONS, DEFLAGRATIONS, AND EXPLOSIONS 

There are important distinctions between the terms detonation, deflagration, and 
explosion. A deflagration is a subsonic combustion wave, behind which the flow is 
accelerated with a significant reduction in density and approximately constant pressure 
(actually a slight reduction). By contrast, a detonation is a supersonic combustion wave, 
behind which the flow is decelerated with higher density and a dramatic rise in pressure. 
[Ref. 7] 

The term explosion refers to an exothermic reaction in which the rate of heat 
release exceeds the rate of heat removal, thereby accelerating the reaction exponentially 
and creating a rapid expansion of the working fluid. The heat generation and pressure 
rise are very rapid, but no combustion wave is required for the process to occur. 
Explosions are thus a "volume reaction” vice a wave phenomenon. [Ref. 5] 

B. THE HUGONIOT CURVE 

The one-dimensional mass, momentum, and energy conservation equations (for 
an adiabatic system in a constant area duct with no body forces) can be combined to yield 
the Hugoniot relation: 

h, ~h -17,-a{2+4] (2) 
Pi P2 
where h is the total (i.e., thermal plus chemical) enthalpy, P is the pressure, and p is the 
density [Ref. 7]. Alternatively, in terms of internal energy, the Hugoniot relation can be 


written as: 


fy SB, Se aa (3) 


Pi, Po 


where e 1s the total internal energy [Ref. 5]. The Hugoniot curve, Figure 3, is a plot of P> 


versus 1/* 5 for a given heat release. 
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Figure 3. Hugoniot Curve from Ref. [7] 

For an adiabatic process, the Hugoniot curve would pass through the origin 
(point A in Figure 3), representing a normal shock (upper branch) or expansion (lower 
branch). With heat release, the curve shifts up and to the right, with greater heat release 
Causing a greater shift. The upper branch, then, generally indicates a supersonic 
combustion process (1.e., detonation) while the lower branch indicates a subsonic heat 
release process (i.e., deflagration). 

For diabatic flow in a constant area duct, the conservation equations also lead to 


the Rayleigh relation: 








Puy = 1 = (4) 


where u is the velocity of the unburned gases relative to the detonation wave [Ref. 7]. In 


terms of Mach number, this becomes: 
Lo (5) 


where y is the ratio of specific heats and M; is the Mach number of the unburned gases 
relative to the detonation wave (or, alternatively, the Mach number of the detonation 
wave propagating into the unburned gases) [Ref. 7]. 

Rayleigh lines passing through the origin of the Hugoniot curve (e.g., the dashed 
lines in Figure 3.) serve to divide it into regions representing different physical 
phenomena. The points of tangency are known as the Chapman-Jouguet points, or C-J 
points. Since a Rayleigh line with a less negative slope could not pass through both the 
detonation branch of the Hugoniot curve and the origin, the upper C-J point (U) yields the 
minimum detonation wave speed [Ref. 7]. It can also be shown that this point 
corresponds to the point of minimum entropy on the Hugoniot curve [Ref. 7]. 
Conversely, the lower C-J point (L) results in the maximum deflagration wave speed and 
the maximum entropy [Ref. 7]. At both C-J points, it can be shown that the burned gases 
travel at sonic velocity (1.e., M2u = M2, = 1) [Ref. 5]. However, in a detonation, the 
burned gases flow in the same direction as the wave, while in a deflagration they flow 


away from it [Ref. 7]. 


The upper C-J point also represents the only steady-state solution in the 
detonation branch of the Hugoniot curve. Region I of the Hugoniot curve (using the 
notation from Figure 3) is called the strong detonation region. Here P2 > Py so uz < wy 
[Ref. 5]. Since M2 < 1, any rarefaction wave, such as could be caused by heat loss, 
turbulence, or friction, would overtake the detonation front and reduce its pressure 
[Ref. 5]. This would shift the detonation down the Hugoniot curve until it reached point 
U, at which point M2 = 1 and rarefaction waves would no longer be able to catch up with 
the detonation front. Thus the wave would reach a steady-state condition. 

Since a detonation is fundamentally a shock that induces chemical energy release, 
the process initially moves upward along the shock Hugoniot from the origin and then 
(with heat release) travels down along a Rayleigh line to the final Hugoniot [Ref. 5]. A 
Rayleigh line that passes through both the origin and Region II must also pass through 
Region I. To travel down a Rayleigh line from Region I to Region II would require 
additional heat release. However, additional shock-induced heat release would simply 
shift the resulting Hugoniot curve up and to the right, with the final state in Region I of 
the new Hugoniot but approaching point U. Thus Region II is not theoretically 
attainable. (However, extremely fast chemical kinetics combined with physical processes 
that do not quite follow the idealized model above could make weak detonation waves 
possible [Ref. 7].) 

Region II, weak deflagration, is simply the combustion phenomena commonly 
observed under most conditions. Region IV, strong deflagration, would require M2 > 1. 


Since heat addition to a constant-area flow cannot accelerate the flow from subsonic to 
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supersonic, solutions in this region are not possible. Finally, in Region V, P> > P; and 
l/* >> 1/* ;. From the Rayleigh relation (Equation 4), u; would be imaginary, so 
solutions here are not physically possible either. 
C. DETONATION WAVE STRUCTURE 

1. ZND Model 

The classic one-dimensional detonation wave model was determined 
independently by Zel'dovich, von Neumann, and Doring. According to this mode], the 
deronaien wave consists of a planar shock wave followed by a chemical reaction region. 
The shock compresses and heats the combustible mixture which, after an induction 
period, reacts. The reaction actually begins immediately but for reactions that follow 
Arrhenius’ Law and have a relatively large activation energy — such as hydrocarbon 
oxidation — the fraction of material reacted changes very little initially but then proceeds 
to completion at a very high rate [Ref. 5]. This chemical reaction then provides the 
energy to sustain the shock wave. The variation of physical properties through the 


detonation wave is shown in Figure 4. 
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Figure 4. Variation of Physical Properties through a Detonation from Ref. [7] 
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A plot of the ZND model on a Hugoniot Curve is given in Figure 5. 
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Figure 5. ZND Model on Hugoniot Curve from Ref. [7] 

The particular path followed (a, b, c, or d, from Figure 5) depends upon the 
reactants being detonated. The progression from a through d represents progressively 
slower chemical kinetics [Ref. 7]. Path a represents the unattainable limit in which 
reaction occurs everywhere through the detonation wave. This would require the reaction 
to begin at the same instant as — and proceed in parallel with — the pressure and 
temperature rise (i.e., exhibiting no causal relationship). Approaching the other extreme 
(path d), the progressively slower chemical kinetics result in progressively larger pressure 
spikes. This spike is short in duration and rapidly disappears as the chemical reaction 
begins. The peak pressure in the detonation wave is referred to as the von Neumann 


spike [Ref. 7]. 





2. Cellular Structure 

In an actual detonation wave, the ZND structure described above is 
hydrodynamically unstable and only observed under transient conditions. [Ref. 7. 8] The 
chemical reactions behind the detonation do not just direct energy toward the detonation 
front. but also transversely. This forms transverse waves, which give detonations a 
complex three-dimensiona] structure. This structure is cellular in nature and appears as a 
fish-scale pattern on smoked foil. Figure 6 provides an example of this pattern and 


Figure 7 illustrates the cellular structure in an axial view of a detonation front. 





Figure 6. Detonation Fish-Scale Pattern from Ref. [5] 
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Figure 7. Axial View of Detonation Front from Ref. [5] 
The fish-scale pattern indicates the path of triple points in the transverse waves. 


There are actually two triple points (A and B in Figure 8), forming a double Mach 


reflection. [Ref. 5] 
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Figure 8. Double Mach Reflection of Transverse Wave from Ref. [5] 
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In Figure 8, the detonation is propagating to the right. Segment JA is the incident 
shock, MA is the Mach stem, and AB is the reflected shock. The hatched lines show the 
reaction front and the dash-dot lines are slip lines. Cell length (L,), defined as the cell 
dimension along the axis of detonation propagation, is determined by the length of 
segment AD. Cell size (A) (also referred to as cell width or cell diameter) is the cell 
dimension perpendicular to the direction of propagation and is governed by the length of 
segment BC. Segment ABCDE forms the transverse wave, propagating downward at 
approximately sonic speed in the burned gases. [Ref. 5] 

Lee (as quoted in Reference 5) describes the cyclic motion as follows: 

Starting at the apex of the cell at A, the detonation shock front is highly 

overdriven, propagating at about 1.6 times the equilibrium Chapman- 

Jouguet detonation velocity. Toward the end of the cell at D, the shock 

has decayed to about 0.6 times the Chapman-Jouguet velocity before it is 

impulsively accelerated back to its highly overdriven state when the 

transverse waves collide to start the cycle again. For the first half of the 

propagation from A to BC, the wave serves as the Mach stem to the 

incident shocks of the adjacent cells. During the second half from BC to 

D, the wave then becomes the incident shock to the Mach stems of the 

neighboring cells. 

The resulting average wave velocity is very close to the equilibrium Chapman-Jouguet 
value [Ref. 5]. 
The chemical reactions are essentially completed within one cell length, and 


A= 0.6L, (6) 


Therefore, as originally suggested by Shchelkin and Troshin (as discussed in Reference 
9), cell size becomes an indication of the induction zone thickness and thus the chemical 
induction time. Accordingly, cell size increases as the sensitivity of the mixture 


decreases. [Ref. 5] 


Unfortunately, cell size itself is not always clearly defined. When the cell size is 
small (1.e., sensitive mixtures). cell structure is usually quite distinct [Ref. 10]. However, 
as mixture sensitivity decreases. multiple instability modes appear in the transverse 
waves [Ref. 10]. These interact nonlinearly to create an irregular detonation structure 
[Ref. 10]. Within this structure, pressures and temperatures can differ by several times. 
causing induction time to vary by two to three orders of magnitude [Ref. 8]. 
Consequently, only the mixture in the first half of the cell 1s directly ignited by the 
detonation shock wave (Region one in Figure 9), and with the hottest points being ignited 
almost without delay [Ref. 8]. Figure 9 identifies by region the different ignition delays 


within a cell of an upward propagating detonation. 





Figure 9. Different Ignition Regions within a Detonation Cell from Ref. [8] 
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Region two requires the additional compression of a transverse wave for ignition, and 
Region three requires still another shock compression [Ref. 8]. Region four, which only 
exists in a portion of the cells, may not react until much later [Ref. 8]. Thus, the effective 
reaction zone is many times larger than that predicted by the ZND model. This complex 
structure makes determination of a dominant mode difficult, thereby limiting the utility of 
empirical laws based on cell size [Ref. 10]. 
D. PROPAGATION OF A DETONATION WAVE 

1. Geometric Diffraction 

Cell size strongly affects the propagation of a detonation through a geometric 
diffraction. When a detonation leaves a constant-area tube and enters an unconfined 
region containing the same mixture, rarefaction waves propagate radially inward. This 
reduces the mixture temperature, thereby increasing chemical induction time. The 
resulting delay may decouple the reaction zone from the shock. The original tube 
diameter must be large enough to allow transverse wave interactions to establish a 
spherical detonation front before the original detonation core is quenched. Otherwise, the 
source of the transverse waves is removed, thus preventing direct re-initiation of the 
detonation and leaving the reaction to propagate as a deflagration. Since the detonation 
core itself propagates at approximately twice sonic velocity in the burned gases, it travels 
a distance of approximately one tube diameter before being quenched by the radially 
propagating rarefaction waves [Ref. 11]. Thus, re-initiation of the detonation must occur 
within this distance [Ref. 9]. Figure 10 shows a successful transition from confined 


planar to unconfined spherical detonation. 
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Figure 10. Detonation Diffraction into Unconfined Volume from Ref. [12] 

The minimum tube diameter required for a detonation to transition from confined 
planar to unconfined spherical propagation in the same mixture is known as the critical 
tube diameter (d,) [Ref. 5, 9]. According to the empirical law of Soloukhin and 
Mitrofanov (as discussed in Reference 9), and supported by observations of other 
researchers [Ref. 5], 

FEY (7) 

The mechanism governing the expansion of a detonation front 1s the generation of 
additional transverse waves, which leads to the creation of very small new cells. 
Transverse waves may be initiated by small random perturbations in the detonation wave 


or by the smallest cells within some of the largest cells in an irregular structure. In the 


latter mechanism, cells in which the induction length exceeds the critical detonation value 
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can actually generate a transverse C-J detonation if ignited by the collision of two 
transverse waves. These transverse detonations (as well as major transverse wave 
collisions) generate convex overdriven detonation fronts in the adjacent unburned 
mixture. This overdriven front initially consists of very fine cells, which then grow to 
form the main structure as the wave returns to the C-J condition. [Ref. 8] 

Less extreme forms of geometric diffraction can facilitate transmission of the 
detonation by providing preferential locations for re-initiation. Murray and Lee [Ref. 11] 
report on a study in which detonations propagated from a rigid tube into a yielding tube 
of the same diameter. By varying the thickness of the yielding wall, they could control 
the strength of the lateral rarefaction wave, thereby allowing re-initiation in less sensitive 
mixtures than would be required for transmission to an unconfined medium. Here the 
yielding wall becomes a source of reignition nuclei for the re-initiation of detonation. 


The process 1s illustrated in Figure 11. 
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Figure 11. Re-initiation in Tube with Yielding Walls from Ref. [11] 

As with propagation into an unconfined region, the detonation core is still 
quenched at a distance approximately one tube diameter from the end of the rigid tube. 
Here, however, an annular detonation forms near the wall behind the detonation core. 
Although the core detonation is required to generate the transverse waves for reignition, 
once detonation near the wall is initiated it will propagate circumferentially and then 
radially until a planar detonation front is reformed, even if the initial core detonation has 
already been quenched. [Ref. 11] 

Murray and Lee [Ref. 11] give four possible causes for the preferential formation 


of reignition nuclei near the wall. Quoting Reference 11: 
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]) it may be induced by reflection of transverse wave remnants from the 
wall; 

2) a thicker layer of shocked but unreacted gas exists near the wall; thus, 
possibly satisfying a minimum volume prerequisite for reignition; 

3) ahigher temperature may exist near the wall due to the presence of "x" 
shaped compression waves behind the slightly curved shock front; 

4) small ripples or other irregularities in the polyethylene surface may 
give rise to localized shock reflections and thus hot spots. 


Another variation is given by Borisov et al [Ref. 13]. By means of a process 
similar to the above, detonation diffraction into cones of less than 60 degrees also causes 
reignition at the boundary. This allows the use of detonation tubes with diameters 
smaller than the critical diameter if a cone is used to assist transmission. Figure 12 


illustrates the arrangement. 





Figure 12. Detonation Diffraction in a Cone from Ref. [13] 
In Figure 12,¢ is the cone angle, segment AA’ is the core detonation, segment 
AC is the lateral oblique shock wave, segment AB shows the radial extent of the core 


detonation at each axial position, and segment AF shows the boundary between the 
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compressed but unburned mixture and the reaction zone. Under critical conditions, 
detonation is reinitiated at the wall just as the core detonation is quenched (i.e., point A 
reaches the axis). At that point, by definition, the total energy in region ACF must be 
sufficient to reinitiate detonation [Ref. 13]. The volume requirements for initiation 
energy, however, are but one of the considerations. For cone angles greater than 60 
degrees, the lateral oblique shock 1s too weak (leading to excessively long ignition 
delays) to support re-initiation of the detonation at the boundary [Ref. 13]. Figure 13 
shows the required diameter for detonation transmission as a function of cone angle. The 


filled circles represent successful re-initiations and the open circles represent failures. 
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Figure 13. Critical Diameter for Transmission through a Cone from Ref. [13] 
2. Compositional Discontinuity 
Desbordes and Lannoy [Ref. 12] examined detonation propagation with a step- 
drop to a leaner (i.e., less reactive) mixture. Since the C-J pressure, detonation velocity, 
etc., are lower in the second mixture (B), a detonation propagating at C-J conditions in 
the initial mixture (A) enters mixture B overdriven. This results in a smaller initial cell 


size in mixture B, thereby reducing the critical tube diameter. Expectedly, as the 
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sensitivity of mixture B 1s reduced, its critical tube diameter increases. However, the 
ratio of its critical diameter to its critical diameter under C-J conditions actually goes 
down. (If the detonation is not overdriven by more than 20%, Equation 7 still applies.) 
This allows transmission from a smaller tube into an unconfined volume. As the 
detonation propagates through B, cell size grows monotonically until C-J conditions are 
reached. [Ref. 12] 

When the detonation in A reaches mixture B, there is a slight delay (which 
increases with C-J cell size in B) in the transmission of the detonation [Ref 12]. In 1963, 
Shchelkin (as discussed in Reference 11) proposed that if the chemical induction time 
increases by an amount on the order of itself, the periodic cellular mechanism will break 
down and the required coupling of chemical and gas dynamic processes will not resume. 
This theory is supported by Murray and Lee observations of cell size approximately 


doubling prior to recovery from near-critical perturbations [Ref. 11]. 


Zo 
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Iii. PDE THERMODYNAMICS 

A. CYCLE ANALYSIS 

Most current chemical propulsion systems are based on constant pressure 
combustion processes. This becomes a standard with which to compare the PDE, which, 
using a C-J detonation, is closer to a constant volume combustion process. The two key 
differences are the amount of heat addition (based on differences in the method of 
combustion) and whether that heat is added along an isobar or a Rayleigh line. | 

The two cycles were compared for a stoichiometric mixture of JP-10 and air (i.e., 
one mole of JP-10 to 14 moles of oxygen and 52.64 moles of nitrogen). As a basis for 
comparison, state zero was taken to be 0.1 atmospheres (representing ambient pressure 
under potential flight conditions). Isentropic compression to 1.0 atmosphere and 300 K 
yielded chamber conditions (under both flight and laboratory conditions) at state one. 
The Thermo-Chemical Equilibrium Program (TEP) [Ref. 14] was used to determine the 
conditions at state two based on the respective cycle. The expansion process from state 
two to state four was taken to be isentropic and the ideal gas law was assumed to hold 
throughout. The parameters from TEP were the frozen equilibrium values and 
recombination of radical species was ignored. 

The constant pressure combustion cycle, or Brayton cycle, is illustrated in 


Figure 14 and the physical properties at each state are given in Table 1. 
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Constant Pressure Combustion 






Pressure (atm) 


Table 1. Constant Pressure Combustion Cycle 


Specific volume at state one (v;) was found using the known inputs and the ideal 
gas law: 


R 
Py=—T 8 
a; (8) 


where P is the pressure, v is the specific volume, R is the universal gas constant 


(8314 ese x), M is the molecular mass (30.4246 Semel at states zero and one), and 
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T is the temperature. Since the compression from zero to one is isentropic, Up was 
calculated using: 

Py = Fy =C (9) 
where y 1s the specific heat ratio (initially 1.3513, from TEP) and C is a constant. 
Temperature at state zero was then found using the ideal gas law (Equation 8). The 
properties at state four were determined in analogous fashion, using the post-combustion 
ratio of specific heats (1.2523) and molecular mass (28.832 ape from TEP. 


The C-J detonation cycle is illustrated in Figure 15 and the physical properties at 


each state are given in Table 2. 
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Figure 15. Chapman-Jouguet Detonation Cycle 
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Table 2. Chapman-Jouguet Detonation Cycle 


The state zero and one conditions are identical to those for the constant pressure 


combustion cycle. State two was determined by TEP and state four was calculated using 


Equations 8 and 9. 


The results above quantify some of the differences between detonation and 


deflagration discussed qualitatively in Chapter II. The two cycles are plotted together in 


Figure 16. 


Pressure (atm) 





Cycle Comparison 


ee SP ~@ 
ae —— Constant Pressure = 
eee oe ——— C-J Detonation | 
it" a 
ee eee 
fi 


| 
i 


Specific Volume (m3/kg) 


Figure 16. Cycle Comparison 
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Figure 16 highlights the dramatic difference between the two cycles. 
Temperature is somewhat higher for the detonation (the rise is actually mitigated by the 
formation of additional radical species) but the real differences are in the pressure and 
specific volume. The detonation pressure at state two is a factor of almost 18.5 greater. 
In addition, while the constant pressure (deflagration) process increased the specific 
volume by eight times, the detonation reduced specific volume almost in half. The much 
larger enclosed area for the C-J detonation indicates a significant advantage in net cycle 
work. 

B. THERMAL EFFICIENCY 

The heat released in a constant pressure combustion process (with negligible 
change in velocity) 1s simply the difference in enthalpy between the states or: 

Gin = hy —h, wis ea (10) 
where qin 1s the chemical energy released, h is the enthalpy, and c, is the specific heat at 
constant pressure. The values at state two were from TEP (cp of 1.4315 “xg x at 2293 K). 
The initial specific heat (1.0397 ais x) was calculated as the mass-weighted average of 
the JP-10 (1.60 ie x) and the air (1.000 eles x) at 300 K [Ref. 14]. The resulting heat 
release was 2970.5 “/¢s. 

Since the compression and expansion were assumed to be isentropic, the net work 
path integral is easily solved: 


Avi = FM 4 p(y ayy) gu Maar 22 
7s) ae 


+ P(v, —v,) (11) 
| A Nee 


Dy q Pav = 


Zo 


where Wc, 1s the net work of the cycle. Using the values in Table 1, result is 1078.1 mele. 


The thermal efficiency for the cycle is the ratio of net work out to heat added, or 
Nos 
Tipe ae (12) 


where Tn is the thermal efficiency. For a constant pressure combustion process under 
these conditions, the result is 36.3 %. 
The C-J detonation process inherently involves flow so more terms are required 


for the energy balance: 


(13) 


where u is the velocity relative to the detonation front. In this case, gin refers to the total 
chemical energy release, which for a detonation includes energy release in the form of 
flow work. Since the reactants are stationary, the relative velocity at state one is simply 
detonation wave velocity (1786.6 "/,). Other state one parameters are the same as above. 
The values from TEP at state two for specific heat and temperature were 1.4532 ie K 
and 2851 K, respectively. The combustion products flow toward the detonation at sonic 
velocity (985.5 ™/,) but the detonation, in this case, is travelling at 1786.6 "/;. The 
relative velocity at state two, then, is the difference between these. The resulting total 
chemical energy release is 2556.2 = 

The net work is calculated the same as above except the C-J detonation cycle 
follows a Rayleigh line from state one to two instead of an isobar. Thus, the first, third, 


and fourth terms in Equation 11 remain the same, or: 
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Canceling like terms yields 


Wann ie steered -v; )-(2 = 1 | =) (15) 


Ve 





Substituting the values from Table 2 into Equations 11 and 15 gives a net cycle work of 
1441.2 i Then, using Equation 12, the thermal efficiency is 56.4 %. 

This impressive potential performance gain generates the interest in PDEs. Even 
with lower heat release, the C-J detonation process is able to perform more net work each 
cycle for the same fuel consumed. The Bese ilty higher thermal efficiency of the C-J 
detonation cycle, compared to that of the Brayton cycle used in current aircraft and 
missile applications, could lead to dramatic reductions in specific fuel consumption. 
Ultimate PDE employment will depend upon the ability to couple this increase in thermal 


efficiency with a reasonable propulsive efficiency in a practical engine. 
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IV. MATERIALS AND EQUIPMENT 

A. JP-10 

JP-10 (CioHj6) 1s a high-density synthetic fuel for turbine driven tactical missiles. 
It is currently employed in military weapons systems, such as the Harpoon, and is 
approved for shipboard use. Although the fuel is not widely studied outside of 
government-funded research, its composition and properties are well known [Ref. 15]. 

Density, kinematic viscosity, and surface tension all affect the atomization of a 
liquid. Since the atomizers were characterized using water, a comparison of these 


properties for water and JP-10 is given in Table 3. 


Table 3. Atomization-Related Properties for JP-10 and Water 












Combustion-related properties of JP-10 are given in Table 4. 


Table 4. Combustion-Related Properties of JP-10 
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B. COMBUSTOR TUBES 

The pre-detonator, air arms, and main combustor are all stainless steel tubes with 
circular cross section. (See the Appendix.) The critical parameters for the pre-detonator 
are its length (long enough for deflagration-to-detonation transition) and its diameter 

greater than the critical tube diameter). At the same time, the volume should be 
minimized to minimize oxygen requirements. At ten inches long, the pre-detonator 
slightly exceeds the 6-8 inches required to generate a detonation wave in a JP-10/oxygen 
aerosol [Ref. 6]. 

As indicated by Equation 7, critical tube diameter is governed by the detonation 
cell size, which itself is a function of the reactants. Unfortunately, a detailed study of 
JP-10 cell size as a function of equivalence ratio (@) and nitrogen dilution (B, where B is 
zero for oxygen and 3.76 for air) does not exist, but efforts at the California Institute of 
Technology are currently underway to characterize these mixtures. Additional efforts at 
Stanford University and the University of California at San Diego are also currently in 
progress to determine the chemical kinetics of JP-10. Until that is completed, however, 
estimates based upon results for other hydrocarbon fuels must be used. Beeson, et al. 
[Ref. 16] report that increasing carbon content appears to slightly increase detonablility 
(i.e., reduce cell size), although (with the exception of methane) the cell sizes are close 
enough so as to be within experimental error. Based on this, the cell size for JP-10 
should be slightly smaller than that for JP-4. Figure 5 of Reference 16 gives JP-4/oxygen 
(B = 0) cell sizes of 3 mm (for = 1) and 1.5 mm (for g = 1.5). TEP calculations indicate 


that the highest detonation velocity for a JP-10/oxygen aerosol occurs near an 
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equivalence ratio of 2.0. Thus, under those conditions, the cell size should be around 
1 mm, resulting in a critical tube diameter of approximately 13 mm. With an interior 
diameter of 1.5625 inches (or 39.7 mm), the pre-detonator should be more than adequate. 

The four air arms have no particular physical or operational size requirements 
beyond providing enough length for good mixing. They have the same interior diameter 
as the pre-detonator but are 13 inches long (to simplify the atomizer piping layout at the 
head end). 

The main combustor has an inner diameter of 5.25 inches (133.35 mm) and a 
variable length based on the number of segments installed. This length should be 
unimportant for initiation since if a detonation cannot be re-established within one pre- 
detonator diameter, the deflagration-to-detonation transition length would be 
prohibitively long [Ref. 5]. Since the pre-detonator exceeds the critical diameter for 
transition to an unconfined volume of the same mixture, there should be no physical limit 
on the maximum main combustor diameter, as long as the mixture uniformly fills that 
volume. In this case, however, the mixture composition changes at the point of 
diffraction. Nevertheless, while not strictly applicable, Equation 7 provides the only 
guidance currently available concerning transition to a spherical detonation front. The 
minimum diameter for stable detonation propagation is known as the limiting tube 
diameter or quench diameter (d ) [Ref. 9]. The applicable relation is 


A 


do =— (16) 
I 


Here, too, knowledge of the cell size is critical. From Figures 5 and 7a of Reference 16 


and Figure 3 of Reference 17, the cell size for a JP-4/air (8 = 3.76) detonation at an 
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equivalence ratio of 1.0 is between about 45 and 60 mm. The cell size for JP-10/air at 
( = 1 may or may not be slightly smaller. Using the worst-case JP-4/air cell size as a 
conservative estimate yields a limiting tube diameter of about 19 mm. This clearly 
should not be a problem. From an operational standpoint, the size of the main combustor 
should be large to maximize the thrust generated. The length, though, will be limited by 
fuel fill time, which 1s driven by the desired cycle frequency and flight conditions. 
C. ATOMIZERS 

The BETE XA-F-PR200 system was selected as the atomizer due to its clean, 
discrete operation during transients and its ability to maintain low Sauter mean diameter 
(SMD or D32) droplets over a wide range of fuel delivery rates. The atomizer is depicted 


in Figure 17. 





Figure 17. Fuel Injection Element from Ref. [18] 

D. VALVES 

All atomization air and oxygen flow was controlled by individual Peter Paul 
Series 20 Model 22 9.5-watt 24-volt DC two-way solenoid valves, each with a 1/8-inch 
orifice (PN 22Z0573NSGM for air and PN 22Z0840NSGM for oxygen). Nitrogen 
pressure to open the spring-pressurized fuel shut-off cylinders in the atomizers was 
applied through Peter Paul Series 20 Model 23 9.5-watt 24-volt DC three-way solenoid 
valves, each with a 3/64-inch inlet onfice and 1/16-inch outlet orifice (PN 


23Z0341NGV). Two lines supplied main air to each of the four air arms, with the flow in 
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each line controlled by Peter Paul Series 70 Model 72 16-watt 24-volt DC two-way 
solenoid valves, each with a 5/32-inch onfice. All solenoids were controlled by Crydom 
6321 optically isolated relays with 100 microsecond response times. 
E. IGNITION SYSTEM 

The ignition system was a Unison constant 1.4 J capacitive discharge system. 
Thirty percent of the energy was deposited in the 40 us spark, for a power of 10.5 kW. 
The spark plug was a Champion (PN CH34168) for aircraft gas turbine ignition. 
F. PRESSURE TRANSDUCERS 

The pressure transducers used were Kistler piezo-electric transducers (PN 603B1) 
with a 2 us rise time. The charge output went to a Kistler charge amplifier (PN 5010B), 
which converted the measured pressure to a 0-5 volt output. A 540 kHz notch filter was 
used on each amplifier to remove the natural frequency of the transducer diaphragm. 
G. CAMERA 

A Pninceton Instruments intensified thermo-electrically cooled charge-coupled 
device (ITE CCD) provided images with a resolution of 576 by 384 pixels. Shutter time 
was set at 500 ns. The camera controller was a Princeton Instruments ST-138 and the 
pulser was a Princeton Instruments PG-200. 
H. SOFTWARE 

Valve and ignition timing were controlled using a Visual Basic 5.0 GUI, which 


was also configured to monitor system pressures. The layout is presented in Figure 18. 
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Figure 18. Facility Control Layout 


V. FLOW CHARACTERIZATION 

A. ATOMIZER CHARACTERIZATION 

1. Particle Distribution 

An Aerometrics Phase Doppler Particle Analyzer (PDPA) system was used to 
determine aerosol particle size distributions. For convenience, water was used in lieu of 
JP-10, since its atomization characteristics have been found to be very similar [Ref. 19]. 
Sauter mean diameters, volumetric mean diameters, particle velocities, and other spatial 
Statistics were obtained by the system. The argon-ion laser sample-volume was placed 
7 cm (2.76 inches) downstream of the tip of the atomizer to obtain an appropriate sample 
with an allowable obscuration. Tests indicated that neither radial nor axial position 
within the spray cone had a measurable effect on droplet size, although not surprisingly 
particle velocity dropped appreciably with distance from the nozzle and the cone 
centerline. Figure 19 1s a picture of the PDPA apparatus layout, showing the argon-ion 


laser, the optical section, and the computer for displaying and recording the data. 
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Figure 19. Phase Doppler Particle Analyzer System Layout. 
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The observed droplet size distributions as a function of atomization pressure are 
shown in Figure 20. The air and fuel pressures were limited to ranges that produced a 


roughly linear response, thereby allowing a linear curve fit of the data. 
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Figure 20. Droplet Sauter Mean Diameter as a Function of Pressure 
Small droplet diameters are critical for the successful detonation of JP-10 ina 
gaseous oxidizer environment. An approach similar to the one by Tulis and Selman 
[Ref. 20] was taken to model droplet heating, evaporation, and oxidation. It was 
concluded from that analysis that droplets of approximately 10 microns in diameter and 
smaller could be heated and vaporized in the necessary time scales required for 
detonation to occur in a JP-10/oxygen mixture at ambient temperature. Figure 20 depicts 


the attainment of Sauter mean diameters below 10 microns for fuel pressures less than 43 
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psig (0.2965 MPa) and air pressures ranging from 80 psig to 90 psig (0.5516 MPa to 
0.6205 MPa). 

Due to the presence of a pre-detonator tube or air arms immediately following the 
atomizers, further flow characterization was done using 5 inch (12.7-cm) and 10 inch 


(25.4-cm) acrylic extensions attached to the end of the atomizer, as depicted in Figure 21. 





Figure 21. Atomizer with Acrylic Extension from Ref. [19] 

These tubes provided an indication of the degree of uniformity, wall adhesion, 
and fuel agglomeration that takes place once the injector is placed inside the pre- 
detonator tube or an air arm. The simplest condition was the injection of a fuel/air 
aerosol mixture into a tube from the head end. The recirculation patterns resulted in fuel- 
rich conditions at the head end that promoted rapid fuel agglomeration. The flow field 
observed for this setup is depicted in Figure 21 (A). The flow field observed for the air 
arm configuration, depicted in Figure 21 (B), enabled a specified fuel-to-oxidizer ratio to 
be loaded with minimal wall agglomeration. Consequently, this configuration allowed a 
wider range of equivalence ratios to be analyzed. Figures 22 and 23 depict the SMDs 
attained at the exit of the 5 and 10 inch acrylic extensions, respectively. In both 


configurations, the minimum SMD attained was approximately 10 microns for a fuel 
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pressure of 40 psi (0.2758 MPa) and an air pressure of 88 psi (0.6067 MPa), with SMD 


ranging from 10-15 microns for all fuel/air pressure ratios. 
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Configuration: Atomizer with 10" pre-det tube 
Location: tube exit 
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Figure 23. Droplet SMD with 10-inch Tube as a Function of Pressure 


2. Flow Rates 

a. Fuel Mass Flow Rate 

Mass flow rate of liquid through the atomizer was determined using a 
pressurized column of water. The nitrogen pressurization system was the same as that 
used for the fuel. Nitrogen pressure was varied from 39 psig to 45 psig in 2-psi 
increments. Atomization oxygen pressure was varied from approximately 80 to 90 psig 
at each nitrogen (1.e., water) pressure. The pressure transducers and the solenoid valve 
control were the same as those used on the pre-detonator. 

For steady-state flow, the atomizer flow was recorded on video, with the 
spray duration taken from the time annotation at the commencement and cessation of the 
flow. Spray duration was typically maintained around 25 seconds. After cessation, the 
change in water column volume was measured using a graduated cylinder. The change in 
volume divided by the duration of the spray was the volumetric flow rate of the liquid. 
Multiplying by the density yielded the mass flow rate. The mass flow rate for each 
oxygen pressure at each fuel pressure was plotted versus oxygen pressure. A linear curve 


fit was applied to the data from each fuel pressure. The results are shown in Figure 24. 
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Figure 24. Atomizer Fuel Mass Flow Rate 

A similar method was used to measure the mass flow rate at different duty 
cycles. Pulses of 500 ms, 250 ms, and 100 ms were examined using the timing in the 
pre-detonator control code. The atomizer was pulsed at the pressure combinations above 
until several seconds of total on-time had been achieved. Then the change in water level 
was measured and the mass flow rate calculated as above. Transient effects were found 
to be significant, with each shorter pulse increment increasing mass flow rate on the order 
of 20% compared to the last. In addition, the oxygen pressure exhibited a roughly 
exponential decay of approximately 3 psi over the first 100 ms of the flow. Since the 
pulses used for detonation are typically a factor of two to four shorter than that, a method 
of estimating mass flow rate from initial fuel and oxygen pressures would have to be 


experimentally determined for each individual pulse duration selected. Instead, it was 


decided to calculate instantaneous mass flow rate every time the atomization oxygen and 
fuel pressures were updated in the code, using values interpolated from the linear curve 
fits of the steady-state data. 

b. Oxygen Mass Flow Rate 

The mass flow rate of oxygen through the atomizer was calculated using 
an evacuated cylinder and the ideal gas law. One end of the cylinder was fastened to the 
nozzle of the atomizer and the other end to a pressure transducer and a vacuum pump. 
The volume of the cylinder and connective tubing (2648 ml) was measured using water. 

After the cylinder was evacuated, the vacuum pump was isolated and 
steady spray was initiated until pressure in the flask reached 15 to 17 psia. Cylinder 
pressure was plotted versus time and two points were selected from the linear range of 
the plot. The change in volume of the water column was measured using a graduated 
cylinder, as above, and both the water and the oxygen were at ambient temperature 
(measured using a thermocouple). Then the ideal gas law was used to calculate oxygen 
mass flow rate: 
R 
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where R is the universal gas constant (8314.3 “/Apmoik) and M o, 18 the molecular mass 


of oxygen (31.9988 “8/comot). 
This was repeated for several oxygen and fuel pressures over the range 
above. The mass flow rate for each oxygen pressure at each fuel pressure was plotted 


versus oxygen pressure. Over this range, oxygen mass flow rate was found to be 
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practically independent of fuel pressure, so a linear curve fit was applied to all the data. 


The results are given in Figure 25. 
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Figure 25. Atomizer Oxygen Mass Flow Rate 

B. EQUIVALENCE RATIO 

1. Pre-Detonator 

Equivalence ratio is defined as the ratio of the actual fuel-to-oxidizer mass ratio to 
the stoichiometric fuel-to-oxidizer mass ratio. TEP provided the stoichiometric oxidizer- 
to-fuel ratio for the combustion of JP-10 in oxygen (3.2882). 

Instantaneous values of the equivalence ratio of the oxygen-and-atomized-fuel 
mixture being injected were calculated by entering the measured oxygen pressure in the 
linear curve fits of oxygen and fuel mass flow rates (Figures 24 and 25), interpolating 


fuel mass flow for the measured fuel pressure, and multiplying by the stoichiometric 
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oxidizer-to-fuel ratio from TEP. Equivalence ratios in the pre-detonator are illustrated in 


Figure 26. 


Equivalence Ratio 
(Pre-D etonator) 


39 psig 
—_ =e» = 41 psig 
43 psig 
—— ——~— — 46 psig 


Ao) 
i 
© 
ao 
= 
2 
© 
= 
S 
a 
Lu 


82 84 
Atom ization 02 prec cure (psig) 





Figure 26. Pre-Detonator Equivalence Ratio 

2. Main Combustor 

The equivalence ratio in the main combustor was essentially the same as that of 
the air arms. Fuel mass flow rate through each atomizer was determined using the same 
plot as for the pre-detonator, except with atomization air pressure as the independent 
variable. Here, however, there are two oxidizer flows that must be considered: 
atomization air and main air. 

The molar flow rate of air through the atomizer is assumed to be the same as that 
for oxygen. Thus the mass flow rate of air through the atomizer was calculated using the 


same data as in Section A, with only the value of the molecular mass changed 
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(M,,, = 28.966 “/Acgmot). The main air mass flow rate in each air arm was determined 


using volumetric flow rate curves supplied by the valve manufacturer. Multiplying the 
volumetric flow rate at a given inlet pressure by the density of air at standard temperature 
and pressure yielded the mass flow rate. This constant was added to the linear curve fit 
of atomization air mass flow rate to yield total air mass flow rate in each air arm versus 
atomizer air pressure. This was performed for main air inlet pressures of 100 psig and 
150 psig. Over the pressure ranges considered here, main air flow is approximately an 
order of magnitude greater than atomization air flow. 

Since the only significant oxidizer in air is oxygen, the stoichiometric oxidizer-to- 
fuel mass ratio 1s the same as that for the pre-detonator. However, since only 20.95% of 
air (by mass) is oxygen, the actual mass flow rate of air must be multiplied by this 
fraction. Instantaneous values of the equivalence ratio of the air-and-atomized-fuel 
mixture entering the air arms were calculated by entering the measured atomization air 
pressure in the linear curve fits of air and fuel mass flow rates, interpolating both fuel 
mass flow for the measured fuel pressure and air mass flow for the measured main air 
inlet pressure, and multiplying by the stoichiometric oxidizer-to-fuel ratio from TEP. 


Equivalence ratios in the main combustor are presented in Figures 27 and 28. 
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Figure 28. Main Combustor Equivalence Ratio (150 psig Main Air) 
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VI. RESULTS 
A. TESTING 

Initially, the pre-detonator was tested to determine optimum conditions. It was 
quickly determined that the presence of the 5-inch diameter main combustor tube had a 
noticeable effect on pre-detonator operation. Operating stand-alone (i.e., exhausting to 
atmosphere), the best pre-detonator operation occurred with a fuel fill time of 22 ms, with 
oxygen flow beginning 3 ms earlier. Ignition was delayed 20 ms after securing 
fuel/oxygen injection to allow mixing. The best operation of the pre-detonator when 
exhausting into the 5-inch diameter tube occurred with fuel fill times of 26-30 ms and 
total oxygen delivery time of 100 ms, where oxygen pressure at the atomizer was 
between 83.5 and 84 psig and fuel pressure was 43 psig. This condition resulted in 
average wave speeds of over 2300 "/; (determined using the time between pressure spikes 
at transducers with a known separation distance). 

The pre-detonator was then tested with main air flow in the air arms. Main air 
was supplied to the solenoid valves at approximately 100 psig and pre-detonator oxygen 
and fuel pressures were maintained as above. Total oxygen delivery time and ignition 
delay were kept at 100 ms and 20 ms, respectively, and main air flow was secured 30 ms 
before ignition. Under these conditions, fuel fill trmes from 30-40 ms yielded average 
wave speeds over 2000 "/;, with detonations occurring in three or four pulses out of five. 
A fuel fill time of 37 ms resulted in four of five pulses detonating and each of the 


detonations being overdriven. 
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A sample timeline illustrating the various timing parameters is presented in 


Figure 29. 
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Figure 29. PDE Test Timeline 
In Figure 29, zero is the start of the cycle. The first three bars (top to bottom) indicate the 
periods of main air flow, atomization air flow, and main fuel flow, respectively. The last 
two bars indicate pre-detonator fuel flow and atomization oxygen flow. To prevent the 
formation of large droplets when initiating fuel flow, atomization air or oxygen flow was 
always established a few milliseconds earlier. The vertical line in the lower nght hand 
corner of Figure 29 shows the point of ignition. 

Finally the entire system was tested with fuel in the main combustor. With the 
pre-detonator and main air conditions as above, fuel fill times through the air arms were 
varied from 98-138 ms with atomization air pressure held at 84 psig. No detonations 
were generated in the main combustor. The injected fuel could be seen outside the end of 


the main combustor, so main fuel fill time was not the issue. 
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It was difficult to determine the source of the problem due to the number of 


variables. Table 5 indicates the ranges over which the key parameters were varied. 


Table 5. Testing Ranges of PDE Parameters 
















Range 


80 - 90 psig 
40 - 43.5 psig 























Although every combination of variables could not be tested, a sufficient range 
was tested to determine that something beyond the scope of these parameters could be 
causing the problem. In addition, several specific steps were taken to address possible 
problems. Delays between flow events and ignition were increased to try to improve 
fuel/oxidizer mixing. Long main combustor air purges were used to remove exhaust 
gases and prevent ignition of the mixture by residual combustion in the tube. Long pre- 
detonator oxygen and fuel/oxygen purges were used to remove exhaust gases from the 


pre-detonator and increase the oxygen concentration at the head end of the main 
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combustor. None of these were effective in enabling the formation of a detonation wave 
in the main combustor. 

Since flowing water through the atomizer produced droplet size distributions 
almost the same as those of JP-10, it was assumed that the volumetric flow rate would be 
the same as well. However, subsequent testing has demonstrated that even with the same 
droplet size distribution, the volumetric flow rate for JP-10 is approximately 20% higher 
than that for water. Thus all tests were conducted at a higher equivalence ratio than 
predicted. However, since equivalence ratio was varied over such a wide range, this did 
not substantially affect the overall results. 

The vapor pressure of the fuel was increased by heating the fuel up to 204°C. 
From Reference 14, JP-10 vapor pressure increases from 0.31 kPa at 20°C to 21.0 kPa at 
120°C. Extrapolating this out to 204°C yields a value of approximately 175 kPa. 
Unfortunately, even this could not generate a detonation in the main combustor. 
However, the increased fuel mass flow rates at these elevated temperatures were not 
accurately determined, so although the equivalence ratio was varied over a considerable 
range so as to examine both lean and rich mixtures, the actual values were not known. 

The installation of a 0.625-inch thick cylindrical quartz section allowed the 
presence of fuel in the head end of the main combustor to be confirmed optically via laser 
transmission. Testing was conducted to determine the pre-detonator fuel fill time (at 
different fuel and atomizer oxygen pressures) required for fuel to reach the main 


combustor. Over our pressure ranges, the longest fuel fill time required was 57 ms. 
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Finally, the Princeton Instruments CCD camera was used to attempt to determine 
what was happening when the detonation wave exited the pre-detonator. The gate width 
was set at 0.5 us and the gating itself was triggered (with a delay) from a pressure 
transducer in the pre-detonator. Figure 30 illustrates a typical detonation failure (i.e., a 


simple deflagration). 
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Figure 30. Failed Detonation 
For comparison, Figure 31 illustrates a detonation from the pre-detonator propagating 


into the main combustor. 
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Figure 31. Detonation Propagating from Pre-Detonator into Main Combustor 
The core detonation is clearly visible, as is the detonation failure around the perimeter 
due to the rarefaction waves propagating radially inward. This shows that the detonation 
front did not transition from planar to spherical. Without re-initiation of detonation at the 
tube wall, the detonation will fail. 
In an attempt to image a successful detonation transmission into the main 
combustor, ethylene and oxygen were used. Although detonation was successfully re- 


initiated, the event shattered the quartz, as shown in Figure 32. 
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Figure 32. Shattered Quartz Section 

Extremely high pressures can be observed at re-initiation, occasionally exceeding 
3000 psi. 
B. FAILURE THEORIES 

Several factors complicate analysis of the detonation propagation in this study. 
Diffraction into an unconfined volume and propagation into a mixture of different 
detonability are complex events individually but having both occur simultaneously 
greatly increases the complexity. Furthermore, aerosol detonations have not been 


thoroughly studied and the effects of droplet size distribution, mixing, evaporation, and 
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two-phase combustion add another level of complexity over gaseous detonations. Lastly, 
detailed chemical kinetics for JP-10 are just now being studied and have not yet been 
analyzed to a level which could aid in PDE design. 

The resulting uncertainties are crucial. Since it is generally accepted that all 
detonation parameters with dimensions of length are governed by cell size, the fact that 
JP-10 cell size as a function of nitrogen dilution and equivalence ratio is not known 
becomes very significant. If the cell size of a JP-10/air aerosol detonation at an 
equivalence ratio of 1.0 were much larger than assumed, the main combustor could be 
smaller than the limiting tube diameter. Below this critical limit, even a single-headed 
spinning detonation would not form [Ref. 9]. It seems highly unlikely that JP-10 cell size 
would be so much larger than that of JP-4, but the example of methane demonstrates that 
there can be extreme differences between hydrocarbon fuels. 

Figure 30 shows a detonation wave leaving the pre-detonator that does not 
transition from planar to spherical when it enters the main combustor. There are several 
possible causes for this including non-uniform fuel distribution at the head end of the 
main combustor, a mixture outside detonability limits, an excessive fraction of the fuel in 
droplets too large to support detonation, JP-10/oxygen cell size greater than 3.05 mm, or 
the empirical law (Equation 7) simply may not hold in this case. 

Another potential problem is the particular geometry at the end of the pre- 
detonator. As indicated in the Appendix, the diameter of the pre-detonator actually 
increases near the exit, from 1.5625 inches to 2.0 inches. This increases the effective 


diameter at the point the detonation enters the main combustor, thereby improving the 
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potential for transition to a spherical wave front or, alternately, allowing greater time to 
re-initiate detonation at the boundary before the core detonation is quenched. However, 
this minor initial diffraction occurs only 0.4375 inches (11.1125 mm) from the plane of 
the main combustor. Murray and Lee [Ref. 11] report that a detonation requires 
approximately 12 cell widths of travel to recover from a boundary-induced disturbance. 
That implies that if the JP-10/oxygen cell size in the pre-detonator exceeds 0.926 mm, the 
initial diffraction actually weakens the detonation entering the main combustor. 

Since the detonation is occurring in a two-phase mixture, the liquid droplets 
(especially the larger ones) almost certainly burn later than the vaporized portion of the 
fuel. This would create a complex and irregular cellular structure. Edwards and Thomas 
(as discussed in Reference 10) suggest that detonations with such structure could require 
a diameter of more than 10-13 cell widths for successful transmission across an abrupt 
expansion. Benedick, et al [Ref. 17], contend that relatively insensitive fuels, such as 
hydrocarbon propulsion fuels, cannot even achieve detonation unless the fuel has 
vaporized. Figure 33 shows the approximate distance behind the detonation front at 


which fuel droplets of varying size vaporize. 
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Figure 33. Fuel Vaporization Behind the Detonation front From Ref. [21] 
For the fuel droplet distributions achieved here, (SMDs of less than 10 microns) the 
results would be close to the infinite acceleration curve. Since the distance behind a 
detonation at which chemical energy release can contribute to the wave is generally taken 
to be approximately 1 cm, a considerable fraction of the fuel is either detonated before 
complete vaporization or simply burns well behind the wave. 

Since the JP-10/oxygen detonation is much more energetic than JP-10/air, 
detonations in the pre-detonator would become overdriven detonation waves in the 
JP-10/air. The more highly overdriven the detonation, the smaller the cell size and the 
more regular the structure [Ref. 10]. This is an important premise upon which this whole 
testing program was based. Table 6 gives an indication of the degree to which the 


detonation entering the JP-10/air aerosol was overdriven. 
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Table 6. Detonation Parameters by Mixture 
Table 6 uses C-J detonation values from TEP. The last row in the table shows detonation 
velocity normalized by C-J detonation velocity in the main combustor (i.e., JP-10/air 
aerosol). Since the pre-detonator was typically operated at an equivalence ratio between 
1.5 and 2.0, the detonation entering the main combustor was significantly overdriven. 
Unfortunately, this was insufficient to generate a spherical detonation front in the JP- 
10/air aerosol. 

The other mechanism for direct initiation of a detonation in the main combustor is 
through the creation of reignition nuclei at the boundary. Several factors were present 
that eliminated this possibility. First and foremost, the ratio of main combustor diameter 
to pre-detonator diameter (see the Appendix) made this unlikely, if not impossible. The 
distance from the pre-detonator to the wall of the main combustor is only slightly less 
than the diameter of the pre-detonator, a ratio of approximately 0.8. For both JP- 
10/oxygen and JP-10/air, the detonation wave propagates at just over 1.8 times sonic 
velocity in the burned gases, thereby travelling 0.9 pre-detonator diameters before being 
quenched. Thus the oblique shock barely reaches the tube wall before the core 


detonation is quenched, leaving almost no time for the formation of reignition nuclei. 
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The only reason the oblique shock has time to reach the wall at all is because the pre- 
detonator tube diameter was expanded from 1.5625 inches to 2.0 inches near the end. 
However, as discussed above, this diffraction occurs too close to the end of the tube to 
allow the detonation to fully recover from the disturbance, so the detonation entering the 
main combustor is somewhat weakened. Furthermore, due to purging of the pre- 
detonator before each cycle, the oxygen concentration in the main combustor is probably 
slightly higher on the axis than at the wall. This would increase sonic velocity in the 
burned gases and thereby accelerate the rarefaction waves quenching the core detonation, 
thus allowing even less time to re-initiate detonation at the boundary. 

Compounding the geometry problems hindering reignition is the initiation delay 
(noted by Desbordes and Lannoy [Ref. 12]) before detonation transmission into the less 
sensitive air mixture. Together these factors would seem to ensure that the core 
detonation would be quenched well before reignition nuclei could be formed at the wall. 
Even without the geometry concerns, there remains the question of whether a detonation 
can be re-initiated at all if the chemical reaction time doubles going from one mixture to 


another. 
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VII. CONCLUSIONS 

A. SUMMARY 

Several complex conditions existed for this testing. A two-phase mixture is by 
definition non-homogeneous, thereby making detonation dependant upon droplet size, 
mixing and vaporization. Moreover, this would tend to produce a multi-modal detonation 
structure, consisting of large and irregular cells. Propagation of such a detonation across 
a geometric diffraction with simultaneous composition change is an extremely 
complicated process. 

Under the conditions tested here, the pre-detonator was shown to be insufficient 
(in energy, diameter, or both) to allow transition from a confined planar detonation to an 
unconfined spherical detonation. At the same time, the main combustor diameter was too 
large to allow re-initiation of detonation at the tube wall. The change in chemical 
kinetics between the mixtures (as well as likely fuel and oxygen non-uniformity) further 
hindered the transmission process. 
B. RECOMMENDATIONS 

Tests should be performed to clarify conditions at the head end of the main 
combustor. First of all, the distributions of both fuel and oxygen at the end of the pre- 
detonator and in the head end of the main combustor should be mapped axially and 
radially. Since JP-10/oxygen detonations are highly repeatable while JP-10/air 
detonations have been as of yet unattainable, detonability of JP-10 as a function of 
oxygen concentration could be tested. Also, an evitiator can be used to vary the degree of 


vaporization of the fuel before detonation. Beginning with fully vaporized JP-10, fuel 
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droplet size should be gradually increased until the maximum Sauter mean diameter that 
will support detonation in a JP-10/air mixture is determined. Since the shape of the 
droplet distribution remains essentially unchanged, a larger SMD directly corresponds to 
a smaller fraction of fuel vaporized. Finally, the ratio of main combustor diameter to pre- 
detonator diameter should be reduced, to allow detonation re-initiation at the wall of the 
main combustor. Once detailed chemical kinetics and cell size information becomes 
available, the relative sizes of the pre-detonator and the main combustor can be optimized 
to minimize oxygen requirements. Alternatively, the method of Desbordes and Lannoy 
[Ref. 12] could be attempted, wherein the detonation propagates into the less detonable 


mixture before undergoing diffraction into an unconfined volume. 
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